THE AUBRY-MATHER THEOREM FOR DRIVEN 
GENERALIZED ELASTIC CHAINS 



SINISA SLIJEPCEVIC 



Abstract. We consider uniformly (DC) or periodically (AC) driven general- 
ized infinite elastic chains (a generalized Frenkel-Kontorova model) with gra- 
dient dynamics. We first show that the union of supports of all the invariant 
measures, denoted by A, projects injectively to a dynamical system on a 2- 
dimcnsional cylinder. We also prove existence of ergodic invariant measures 
supported on a set of rotationaly ordered configurations with an arbitrary 
(rational or irrational) rotation number. This shows that the Aubry-Mather 
structure of ground states persists if an arbitrary AC or DC force is applied. 
The set A attracts almost surely (in probability) configurations with bounded 
spacing. In the DC case, A consists entirely of equilibria and uniformly slid- 
ing solutions. The key tool is a new weak Lyapunov function on the space of 
translationally invariant probability measures on the state space, which counts 
intersections. 



1. Introduction 



(1.1) - Uj (t) = - V 2 («,-! - V 1 {u j {t),u j+1 {t)) +F(t), j G Z 



We consider the equation 
d 

—i 
dt 



such that V : M 2 — > R is C 2 , periodic (V(u + 1, v + 1) = V(u, v)), and satisfies the 
twist condition Vi2(tt,u) < — 5 < for some fixed S G K. 

We assume that F is either constant (the DC case) or 1-periodic in t (the AC 
case). In the AC case, we further assume that V, F are real analytic. 

1.1. The background and motivation. The case F — is the standard setting 
of the well-known Aubry-Mather theory ([3], jTB] , Section 9.3), as the solutions of 
du/dt — correspond to orbits of an area-preserving twist map fv on a cylinder 
whose generating function is V. The Aubry-Mather theory describes minimizing 
configurations/orbits and measures, and so simultaneously characterizes physically 
important "ground states" of (|1.1[) and describes "KAM-circles" and what is left 
of them for fv- 

Our goal is to contribute towards the description of the dynamics of (jl.ip on the 
entire state-space, allowing a fairly general family of initial conditions. We are, how- 
ever, inspired by results on asymptotics of one-dimensional reaction-diffusion equa- 
tions ([S], [H]), their braid dynamics ([S]); stability of synchronization phenom- 
ena ([IH]); m the Hamiltonian dynamics setting the Mather's measure-theoretical 
approach to minimizing orbits and its extensions ([H]); PDE extensions of Aubry- 
Mather theory ([H]); results on spatio-temporal entropy of extended PDE's ([2D]. 
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[35]); all related to what is actually known on the behavior of (jl.ll) either numeri- 
cally, physically ([7]) or analytically ([1], [5], [H], [22], [23]). Before stating precisely 
the actual results, we will try to put them in a wider context of these topics. 

In the "standard" case V(u, v) = (u - v) 2 /2 - W(u), where W(u + 1) = W{u) 
and W : ffi — >• R is C 2 , the equation is a spatial discretization of the scalar 

reaction-diffusion equation 

u t = u xx + W (u), 
and it is closely related to a more general reaction-diffusion equation 
(1.2) u t = u xx + f(x,u,u x ), 

where / is 1-periodic in the first coordinate. Both (jl.ip and (ll.2p are scalar and 
order-preserving (or monotone, or cooperative, see [27]). The well-known Fiedler 
and Mallet-Paret Poincare-Bendixson theorem [5] showed that the asymptotics of 
(|I.2j) on any bounded domain is 2-dimensional, which is a basis of a relatively 
complete understanding of dynamics of (|1.2[) on bounded domains f |15j and refer- 
ences therein). The key tool they applied is the intersection-counting function 
(u,v) — > z{u — v), which is a discrete Lyapunov function. That means that 
t h-> z(u(t) — v(t)) is non-increasing along two solutions of (II.2j) . Furthermore, 
if u(to) — v(to) has a singular zero (we also say that u(to) and v(to) intersect non- 
transversally), z(u(t) — v(t)) strictly drops at t = to- 

These ideas have not in our knowledge been applied to the equation (|1.2p on 
the entire K, considered as an "extended system" (i.e. without assuming u(x) — > 
as x —¥ oo). The reason is that the intersection-counting function z is typically 
not finite when it, v are defined on the entire K, thus all bounded-domain tools 
fail. We believe that the right approach in the extended-domain case is to consider 
the extension of the intersection-counting function z to the space of translationally 
invariant measures on the state space. 

In this paper we pursue this idea for a technically simpler discrete-space system 
(jl.ip , and consider a function 

(a*\m 2 ) 1 — > z(m\m 2 ) 

which counts the average number of intersections of configurations with respect to 
two translationally invariant probability measures on the space of configurations 
{ui)iez- We will show that t h-> Z (n 1 (t) , /j, 2 (t)) is a weak Lyapunov function (for 
details see Section 0, where //(£), A* 2 (i) are evolutions of measures with respect to 
(jl.ip . We believe these ideas can also be applied also to (|1.2j) and give fuller insight 
on the dynamics of (|1.2p on unbounded domain. In particular, we conjecture that 
all our results (with appropriate rephrasing) hold also for p.2j) . 

Van der Berg, Ghrist, Vandervorst and Wojcik studied in several papers topo- 
logical " braid" properties of intersections of solutions of (|1.2p on bounded domains 
and its spatial discretization ( |13j and references therein). Our approach may lead 
to extension of their ideas to unbounded domains. 

It has been noted that synchronized solutions of the extended reaction-diffusion 
equation in more spatial dimensions (> 2) are often stable and locally attracting 
(see e.g. [10] for recent results and references). We argue that it is possible to study 
this phenomenon also in the spatial dimension 1. As compared to we also allow 
unbounded solutions Ui{t) (by identifying u and u + 1 and considering the dynamics 
in an appropriate quotient space), and also study non-autonomous systems (AC 
forcing). We eventually reduce our dynamics to a 2d flow (DC case) or map (AC 
case) with non-trivial dynamics (the spatially Id dynamics in [lOj reduces to a 
trivial Id flow). A heuristic explanation for stability of Id synchronization is now 
as follows: the synchronized solutions are the solutions in the support of invariant 
measures \x such that the self-intersecting function Z([i,fj,) is equal to zero (we 
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consider measures invariant both for spatial translations and time evolution). As 
the function t i— > Z (fi(t) , /i(t)) is non-decreasing, such solutions are expected to be 
at least Lyapunov stable. We develop a rigorous approach to this idea. 

The Aubry-Mather theory of ground states of in the case F — has been 
extended to more dimensional lattices and PDEs ([IB] and references therein). In 
our knowledge, generalizations typically deal with existence of ground states with an 
arbitrary generalized rotation number (rational or irrational). In the PDE setting, 
this leads to construction of various solutions of stationary PDEs. In our view, 
it is somewhat unexpected that the Aubry-Mather stationary structure seems to 
be preserved in a not entirely obvious way for an evolution equation (jl.ip with 
arbitrarily strong, and even time-periodic forcing. 

Finally, let us come back to the actual dynamics of The Frenkel-Kontorova 

model is one of the most studied physical models, as it reasonably realistically mod- 
els a variety of physical phenomena, and is perhaps a simplest model with suffi- 
ciently rich structure allowing lots of "metastable states", various phase transitions 
such as "pinning-depinning" , etc. ([7] and references therein). It has been long 
noted in the physics literature that the dynamics of (11.11) seems to be attracted 
to low-dimensional invariant manifolds, or to spatially synchronize. Baesens and 
MacKay proved that in the DC case spatially periodic initial conditions converge 
either to equilibria or to uniformly sliding solutions ([4], [5]; see below for details). 
Hu, Bambi and Qin extended some of those results to show existence and properties 
of spatially periodic, ordered solutions in the AC case ([H]), and Qin studied exis- 
tence and properties of uniformly sliding solutions with both rational and irrational 
rotation number in the DC case ([22], [23] ). 

Here we attempt to give a broader perspective on the dynamics of (|1.1|) on a 
large state-space, and generalize most of the mentioned results. Firstly, we consider 
dynamics on the space of configurations of bounded spacing (configurations (tii)igz 
such that sup igZ |iti+i — Ui\ < oo), which is more general than the configurations 
of bounded width (i.e. such that sup igZ \vn — ip\ < oo for some p £ R) rigorously 
studied so far. This is important, as it includes configurations connecting two 
configurations with different rotation number, and various other relevant situations. 
We then show existence of a 2d invariant set of configurations in both AC and DC 
case which in a physically meaningful way attracts configurations with bounded 
spacing. Finally, we show existence of ordered configurations with an arbitrary 
(rational or irrational) rotation number. 

1.2. The statement of results. The state space for the dynamics of (jl.ip of 
configurations of bounded spacing will be denoted by X, and defined precisely as 
follows: 

K n = [iid R z , such that sup |iti+i — ui\ < n}, 

X = Kn, 

n=l 

with the induced product topology. We will use the quotient sets X = X/R, 
JC n = K n /R, where Ru = u + 1 is a translation commuting with (II. ip (for details 
see Section [5]). 

The system generates a continuous global semiflow and a local flow on X 
(see Section 4), which we denote by ip (or tp for a given t). Of key importance will 
be the translation T : X — > X defined naturally as (Tu)i = tii—i, commuting with 
(HU) and <p. 

The Borel probability measures on X will be called measures with bounded spac- 
ing. We will in particular study measures with bounded spacing, invariant with 
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respect to both the translation T and the semilflow p> 1 , called (ip, T)-invariant. In 
the DC case, the measures will be invariant with respect to cp* for all t € K, and in 
the AC case with respect to t € Z. Denote the sets of such measures on JC n , X by 
£(7C„), Let 

(1.3) .4 = U Me£(A .)Supp(/i), 

where supp(/i) denotes the support of /i. Our first main result is the following: 
Theorem 1. There is a continuous infective projection ir : A — > T 1 x R. 
The injective projection will be given by 

(1.4) tt(u) = (u ,ui - u a ). 

As the projection ir is injective, we can pull the translation T, and the semiflow 
ip f ' in the DC case, or the map ip 1 to the 2-dimensional set 7r(X). We introduce the 
following notation: 

X = n(X), 

/ 1>5 \ ^„ = tt(/C„), 

T = 7ToTo 7T _1 , 

~i t -1 

</? = TT O ip L O IT . 

Theorem 2. For each n G N, t/ie map T is a homeomorphism of K, n , in the 
DC case pf is a continuous flow on IC n , in the AC case if is a homeomorphism on 
IC n , and T and ip (respectively Tp 1 ) commute. 

We propose to call the pair (T,tp) in the DC case the characteristic flow-map, 
and the pair (T,ip ) in the AC case the characteristic maps. 

We give a detailed example in Section [2] for the case F = 0, and relate our main 
results to the Aubry-Mather theory. In that case in particular, X is the entire 
cylinder T 1 x IR if V is close to integrable, lp is the identity, and T is diffeomorphic 
to an area-preserving twist map whose generating function is V ([16]). 

We do not know whether in general the maps T and </3* can be naturally extended 
to the entire cylinder T 1 x M. Also unlike in the equilibrium case F = 0, we do not 
a-priori know T, Tp . We give some conjectures in Section|51 We can, however, show 
that X C T 1 x IR is a" rich" set for all F , as it contains analogues to Aubry-Mather 
sets of minimizing configurations/orbits with an arbitrary rotation number (called 
also mean spacing) found in the case F = 0, and to Mather's minimizing measures 
([19]). (We postpone definitions of rotation number, a (ip, T)-ergodic measure, a 
rotationally ordered configuration, and totally ordered orbit to Sections [SJ [5]) 

Theorem 3. Given any p 6 K, there exists a (ip,T)-ergodic measure on X sup- 
ported on a set of rotationally ordered configurations with the rotation number p. 

In the DC case, we can completely describe the set of all (<p, T)-invariant and 
ergodic measures: 

Theorem 4. Assume F is DC. Then every ((p,T)- ergodic measure on X is sup- 
ported either on a set of equilibria, or on a single totally ordered periodic orbit. 

The totally ordered periodic orbits in Theorem U are in the literature known 
as uniformly sliding states ([3], [S], [22], [23]). We show in Section that the 
known results on uniformly sliding states follow as a direct corollary from our 
considerations. 

Finally, we address asymptotics of the dynamics (|1.1[) relative to the set A defined 
by (USD. 



AUBRY-MATHER. THEOREM FOR ELASTIC CHAINS 



5 



1.3. The spatio-temporal attractor. We propose a notion of attractiveness 
which is weaker than the usual notions of attractor or union of w-limit sets, but 
for extended systems such as ([l.ip often in practice (physical systems or numerical 
simulations) indistinguishable from them. Imagine a physical situation in which wc 
observe dynamics of on Z (or a very large subset of Z), and that no position 
in space is preferred. Consider also that an initial condition it(0) emerges as a real- 
ization of some spatially homogenous random process. For example, we can choose 
u(0) as a double-infinite sequence so that Ui (or u^+i — Ui) is chosen randomly 
and independently in each step i £ Z from the same bounded measurable subset 
of K. We can study that as analyzing the dynamics (jl.ll) with respect to fi-a,.e 
initial condition, where /x is a T-invariant probability measure on X. We propose a 
notion of spatio-temporal attractiveness, and say that a set A is a spatio-temporal 
attractor, if given any T-invariant probability measure /i, then for /i-a.e. initial 
condition u and with density of times t equal to 1, </?*(«) is arbitrarily close to A. 
More precisely, given any S > 0, we define the Borel probability measure -P^,s on 
X x [0,S] as the product measure fi x As, where As = (l/S , )A|[ 0i s] is the normed 
Lebesgue measure A restricted to [0, S], and define the attractor as follows: 

Definition 1. A set A is a spatio-temporal attractor for on X, if given any 

T-invariant probability measure \i on X and any open neighbourhood U of A, 

(1.6) Urn P^ s ({(u,t), p*(u) G U}) = 1. 

S— >oo 

We finally in Section [TU] we show that the set A described in Theorems QJH is 
indeed a spatio-temporal attractor. 

Theorem 5. The set A defined with \1.3)) is the spatio-temporal attractor for M.l)) 
on X . 

Lastly, we comment on the spatio-temporal entropy of (11.1[) . A corollary of The- 
orem 2] is that the spatio-temporal entropy of (|1.1[) on X in the DC case is zero 
(see e.g. [10], for definitions). A proof follows from the variational princi- 
ple for metric and topological spatio-temporal entropy and our description of all 
spatio-temporal invariant measures; we intend to report details separately. This 
complements well the existing results for PDE systems such as fj 1 . 2|) . Specifically, 
Zelik [52] showed that extended gradient systems (the case F = of (jl.ll) ) have 
spatio-temporal entropy zero, Mielke and Zelik 20 constructed Id non-autonomous 
PDE examples with positive entropy, and Turaev and Zelik [28 proved positive en- 
tropy for a Id complex system. 

1.4. The structure of the paper. In Section[2]we as an example summarize the 
known facts on the dynamics (11.11) in the "extended gradient" case F = 0, and relate 
them to our results. In Sections^ 2]we develop tools on the dynamics of the zero-set 
of a non- autonomous system of infinitely many linear ODEs, and as a corollary on 
dynamics of intersections of a cooperative system of infinitely many ODEs. Here we 
extend some of the known results ( [5] , [25] , [2JJ ) to systems of infinitely many ODEs, 
less smooth then e.g. in [26], and on a general space of solutions (the Frechet space 
of configurations of sub-exponential growth at infinity) . As such, our results could 
be of independent interest. We postpone most of the technicalities to Appendices 
1,2, dealing in particular with some exceptional, "degenerate" situations due to the 
infinite number of equations. In Section [5] we introduce the semiflow on the space 
of translationally invariant measures on X, and show that the average-number-of- 
intersections-counting function Z(/x 1 ,/x 2 ) is non- increasing and strictly decreasing 
if u,v in the supports of /x 1 ,/i 2 have a non-transversal intersection (or, using the 
alternative terminology of the paper, u — v has a singular zero). Once all the tools 
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are in place, the proofs of Theorems [JJ[5] are relatively straightforward and given in 
Sections ICTTU1 including a more detailed discussion of uniformly sliding states. 

2. Example: dynamics without driving and Aubry-Mather theory 

In this section we relate our results to the known results in the case of no driving 
F = 0. Our assumptions on the function V are the standard setting of the Aubry- 
Mather theory. Let us first consider the equilibria du/dt — of (jl.ip . denoted by 
£ C M. z /R, where we again identify u and Ru = u + 1 and consider the quotient 
set. It is straightforward to deduce from the conditions on V that the projection 
7r : £ — >T l xl is a homeomorphism (again considering the induced product topology 
on £). Also the map T : T 1 x K — > T 1 x R is C 2 , diffeomorphic to the twist map 
generated by V ([16], Section 9.3). The diffeomorphism is given by the discrete 
Legendre transform, which is for example identity in the "standard" case 

(2.1) V{u,v) = {v -u) 2 /2 + W(u), 

where W(u + 1) = W(u) and W is a C 2 function. 

It is however not a-priori clear that all the (ip, T)-invariant measures of (|1.1[) 
with F = are supported on £ . The system (|1.1[) is an example of an extended 
gradient system, i.e. a system which has a formal Lyapunov function (in our case 
the formal sum X^-oo V( u i> u i+i))> an d is a gradient system when reduced to a 
finite domain. It has been shown that extended gradient systems can have more 
complex dynamics than what is by LaSalle principle possible for gradient systems 
f [11). [T2l . [25]). including possible existence of non-stationary orbits in cj-limit sets. 

In |25j we however showed the following: 

Theorem 6. All (ip,T) -invariant measures of U.l\) on X in the case F = are 
supported on £. 

The key tool is an observation that the function /i t— > J V (uq, ui)dfj, is a Lyapunov 
function on the space of T-invariant probability measures X ', so LaSalle principle 
must hold for the evolution of T-invariant measures. 

We conclude that Theorems [T][3] do hold in the case F — 0, where the measures 
in Theorem [3] are for example Mather's minimizing measures ([B], |19|). 

Let us briefly discuss the question whether and when it : A -^T 1 x R is a home- 
omorphism, where A is the support of all (ip, T)-invariant measures of bounded 
spacing as defined in (|1.3[) . One can see that tt is a homeomorphism if and only 
if £ C X. This is true because T is diffeomorphic to an area-preserving map, 
so we can pull back a Lebesgue measure and construct (ip, T)-invariant measures 
supported on any closed subset of £ n X. In general, £ C X is true if and only 
if the area-preserving twist map fv generated by V is close enough to integrable. 
For example, in the standard case (|2.1[) . £ C X holds if and only if fv has homo- 
topically non-trivial invariant circles, otherwise one can construct "accelerating" 
invariant measures of fv', i.e. there are T-invariant measures in £ which do not 
have bounded spacing (pHj). 

Finally, let us consider asymptotics. The following slightly stronger version of 
Theorem \5\ in the case F = was shown in and extended in [TT], [T^] for much 
more general systems: 

Theorem 7. Consider the dynamics with F — on X . 

(i) If /J, is a T-invariant Borel probability measure on X , then the uj-limit set for 
fi-a.e. u € X consists of equilibria. 

(ii) Let u £ X and let U be any open neighborhood of £. Then the set of times 
t £ [0, oo) such that u(t) € U has Banach density 1. 
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We can rephrase Theorem [71 and say that the oj- limit set consists of equilibria 
either almost surely with respect to space, or with respect to time. One can not 
omit the " almost surely" part of the statement and still obtain a 2d union of limit 
sets either in Theorem [5] or [JJ as it is most likely infinite-dimensional, analogously 
to the continuous space case ([2T]). 

3. The dynamics of a zero-set 

Here we first consider in general terms properties of the zero-set of a solution of 
a linear, cooperative non-autonomous system of infinitely many ordinary differen- 
tial equations. We show that the number of zeroes in a bounded domain can not 
increase and is strictly decreasing at singular zeroes. We then extend these ideas 
to systems of infinitely many equations, and establish a "zero-balance equation". 
Although some of the ideas presented here are essentially for finite dimensional 
systems already contained in [33], [37], our presentation follows more closely An- 
genent, Fiedler and Mallet-Paret [I], [2], [9] who developed it in the PDE setting. 
In particular, we use the PDE-terms "regular" and "singular zeroes", rather than 
"transversal" and "non-transversal intersections" of v — 0, as we will in a way 
consider the dynamics of the zero-set. 

Let us consider a system of equations 

(3.1) — Wj (t) = aj (t)wj-x + bj {t)w j+ i + Cj (t) Wj , j eZ 

We assume that for some To < Ti, the following holds: 

(3.2) 



aj,bj,Cj continuous, uniformly bounded in j,t <E [Tq,Ti] 



aj, bj > S > for all j, i G [T , Ti] 

The condition a,j,bj > 5 > means that the system (|3.1[) is cooperative f|26j. 
|27j). We can consider (|3.1[) to be a discrete analogue of parabolic differential 
equations considered by Angenent (in real analytic case in [T], in C 2 case in [2]). 
We consider (|3.ip on the space Y of configurations of sub-exponential growth. We 
define it as the set Y C K z of w — (wi)i e z such that 

1- lo SK-| n 

hm sup — = 0. 

j-s-oo \J\ 

Indeed, Y is a Frechet space generated by the family of norms 
(3.3) Mkoo = supexp(-|j|/n)|u|j, 

n a positive integer. 

We first give definitions. We say that (j,t) C Z x [To,Ti] is a zero of a w : 
[To, Ti] — > Y", if the graph of i M> Wi(t) intersects y = in the interval [j, j + 1) (i.e. 

Wj{t) + (w j+ i(t) - Wj(t))x = 

has a solution for x € [0, 1)). We say that (j, t) is a singular zero, if Wj(t) — 0, and 
either of the following three cases holds: 

w j+ i(t) = 0, or 
Wj-i(t) = 0, or 
Wj +1 (t)wj-i(t) > 

(these are the discrete space analogues of w x (t) = 0); otherwise (j,t) is a regular 
zero. We denote by Z, S, R C Z x [T ,Ti] the sets of zeroes, singular zeroes and 
regular zeroes of a given w(t). A singular zero has a degree k, if k is the maximal 
number of consecutive j such that (j, t) is a zero. Assume now (j+1, t), {j+k, t) £ 
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S is a regular zero of a finite degree k, (j, t), (j + k + 1, t) ^ Z. We distinguish two 
types of singular zeroes: 



Note that for singular zeroes of Type I, the degree is k > 2, and for Type II, k > 1 
(see Figure 1 in the Appendix 1). A zero can also be of degree oo, that is w(t) G Y 
such that for all i > io, itfj(t) — 0, or such that for all i < io, tVi{t) — 0, but w 
not identically 0. When discussing zeroes, we will need to consider this, intuitively 
degenerate and quite unlikely case separately with some care. 

Given w £ Y, we define zj(w) to be 1 if w has a zero at j, otherwise 0. Let 
z m ,n ■ Y — > NU{0} be the zero-counting function 



where m < n are integers. 

We say that a solution of (|3 . 1 [) is C 1 if it is C 1 in each coordinate. We can 
now state the discrete analogue of Angenent [2], Theorem D (i.e. for free boundary 
conditions not intersecting zero): 

Lemma 1. Assume w : [Tb,Ti] — » Y is a C 1 solution of \3.1}) satisfying i3.2\) . 
m <n integers. Assume also Wj(t) ^ for (j,t) € {m,n} x [Tq,Ti]. Then 

(i) The number of zeroes z m>n {w{ij) is non-increasing. 

(ii) If w(to) has a singular zero between m and n, then z m>n {w{t)) is strictly 
decreasing at t = to. 

(Hi) If w(to) has a singular zero between m and n, then there exists 5$ > so 
that for all < \S\ < w(to + S) has no singular zeroes between m and n. 

We postpone the proof (as well as all the other proofs in this section) to Appendix 
1 as it is rather technical, and the approach is not related to the rest of the paper. 

We would now like to make general statements on values of z„- hn for arbitrary 
time and space intervals, including those when a zero is "crossing" a boundary 
{m,n} of a segment, and the assumptions of Theorem [T] do not hold. 

We will associate to a solution w : [Tq,Ti] — > Y of (|3.1[) a family of functions 
Ci(w;s,t), i G Z, < s < t, counting the number of times a zero enters minus the 
number of times it leaves the segment i > during the time interval (s, t] (a precise 
definition will be given below). In addition, we define di(w;s, t), i 6 Z, < s <t 
as the number of zeroes "disappearing" at the position i in the time interval (s,i\. 
We will omit the argument w when it is clear from the context. 

As the definitions of Ci , rfj are intuitively straightforward but technically some- 
what subtle, we give it axiomatically, and then prove their existence in Appendix 
1. 

(Al): The zero-balance equation: for all integers m < n and all Tq < s < 

t < T\ the following holds: 



(A2): Translation invariance: Cj,dj are invariant for the translation T, 
that is 



Type I, if Wj(t)wj +k +i(t) < 0, 
Type II, if Wj(t)wj +k +i(t) > 0. 




n-l 



(3.4) 




Cj(Tw; s, t) 
dj(Tw; s, t) 



Cj+i(w;s,t), 
d j+1 (w;s,t). 
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(A3): Regularity of c: If Wj(t) ^ for each t E [to,ti], then Cj(to,tj.) = 0. 
(A4): Cardinality of d: If {(m, to), (m + k — 1, to)} is a zero of degree k, 
then there is So > so that for all < 5 < So, 

E'° 1 r , , f fc — 1 zero of Type I, 

dj(to-S,t + S) = i , J * ' 

J v ' I fe zero of Type II. 

If there are no singular zeros at j for t E [to,*i], then dj(t ,ti) = 0. 
We again postpone the proof of the following to the Appendix: 

Lemma 2. We can associate to each C 1 solution w : [Tq,Ti] —> Y of h3. 1\) satisfy- 
ing \S.2\) a family of functions Cj(w;s,t) E Z and dj(w;s,t) E {0,1,..} satisfying 
(Al)-(Mh 

Now recall the partial order on !L Z : 

w > v if for all j, > Vj. 

The following well-known monotonicity property ([5], can be understood as a 
special case of the discussion above; for the sake of completeness a short proof is 
also in Appendix 1. 

Lemma 3. Assume w : \To,T{\ — >• Y is a solution of \3.1)) satisfying Then 
ifw(0) > 0, then for all t > 0, w(t) > 0. 

4. The dynamics of intersections 

In this section we consider the dynamics of intersections w — u 1 — u 2 of two 
solutions u 1 ,u 2 of (jl.ll) and of derivatives w = du/dt of a solution u of We 
first show that the results of the previous section apply to u 1 — u 2 in both the AC 
and DC case, and to du/dt in the DC case. We then establish local and global 
existence of a continuous semiflow generated by (jl.l[) . Finally in a series of four 
lemmas, we establish the key tool: that the singular zeroes of u 1 —u 2 and du/dt are 
persist for small perturbations in the product topology. Here we treat the cases of 
singular zeroes of finite and infinite degree separately, as the later case is technically 
more different. 

4.1. Existence of semiflow. Recall the definition of K n as the set of all u E IR 2 
such that \uj + \ — Uj\ < n, and fix n E N. 

Proposition 1. Assume w^u 2 ,?! : [Tq,Ti] — > K n are solutions of hl.l}) . Then 



w = u 2 — it 1 in 
some ai,bi,Ci satisfying \3.2\ 

Proof. If w 



all cases, and w = du/dt in the DC case, are solutions of \3. 1\) with 



V 2 (u 2 _ v uj) -V 2 (u}_ 1)U }) 



1 2 5 

U i+l ~ u i+l 

V 2 {u 2 _ 1 ,u 2 ) - V 2 {u 2 _ 1 ,u\) + V^ul u 2 +1 ) - ViKSuLi) 



(naturally extended to ai — — V\ 2 , bi — —V\ 2 , Ci — —V\\ — V 2 



22 



u 2 _ 1 , u\ +1 = uf +1 , u\ = uf respectively). Then by the twist condition and the mean 
value theorem, <Zj, 6, > 8 > 0, and by periodicity of V and a simple compactness 
argument, at{t), h(t), Ci(t) are bounded, thus Q3.2p holds. As u 1 (i),M 2 (i) E K n , the 
difference j M> \uj(t) — u 2 (t)\ grows at most linearly in \j \ , thus M 1 (i) — u 2 (t) E Y. 
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Now let w = du/dt in the DC case. Then by differentiating (jl.lj) we see that 
([XT]) . (EJ2) hold with 

Oi = -Vl2(Ui-l,Ui), 

h = -Vi2(ui,u i+ i), 

c% = -V 2 2(ui-i,Ui) - Vu(ui,u i+1 ). 

a 

We now establish that the equation (jl.lj) generates a smooth semiffow on JC n = 
K n /R, denoted by Lp f and tp. 

Lemma 4. The equation (1.1]) generates a continuous semiflow on ]C n , and a 
continuous local flow on X . 

For each i € Z, the solution t i— > Ui(t) is C 2 ; and if F,V are real analytic, 
t i y Ui(t) is also real analytic, all on all open sets of all t for which the solution 
u{t) G X. 

Proof. It is easy to show by applying the well-known results on existence of solutions 
on Banach spaces, that for u(0) € K. n the solution of (jl.ip exists locally, and 
t M> u(t) — u(Q) is continuous in ^(Z). Thus (jl.ip generates a local flow on X. 

We can show that the solution depends continuously on initial conditions in the 
product topology, by for example establishing local existence of solutions in any 
of the norms ||.|| n ,oo introduced in the previous section (note that these norms 
induce the product topology on /C„ and X). Finally, given a positive integer n, the 
property u £ K, n can be written as 

Tu + n > u > Tu — n, 

which is by Lemma EH and Proposition [T] applied to w(t) = (Tu + n)(t) — u(t), 
w(t) = u(t) — (Tu — n)(t) invariant forward in time. We conclude that if u(0) G /C n , 
then u(t) e K, n for all t > 0. 

Regularity follows from (jl.ip . In the real analytic case we use real analyticity of 
F, V and check inductively that the derivatives satisfy the characterization of real 
analyticity ([H], Proposition 1.2.12). □ 

We will often fix an integer n and write K, omitting its superscript. 

Generally the results to follow regarding the difference u 1 — u 2 will hold in both 
AC and DC cases, and the results regarding the derivative du/dt in the DC case 
only. 

4.2. The AC and DC cases. We first show that singular zeroes of a finite degree 
persist for small perturbations in the product topology. Without loss of generality 
we consider zeroes at to = 0. 

Lemma 5. Assume u 1 (0),u 2 (0) £ JC such that m 1 (0) — w 2 (0) has a singular zero 
of a finite degree k at (m + l,to) 5 {m + k,to). Then for all sufficiently small 
e > there exist open neighborhoods U 1 and U 2 o/u 1 (0), u 2 (0), such that for all 
w 1 (Q),u 2 (Q) e U 1 x U 2 , 

m-\-k 

(4.1) J2 d^-v 2 - -£,+£)> -. 

j=m+l 

Proof. By continuity and Lemma [I] (iii) there exists 6q > so that u x (t) — u 2 (t) 
is not zero on {m,m + k + 1} x [— Sq,5o] and so that u x (i) — u 2 (t) has no other 



AUBRY-MATHER. THEOREM FOR ELASTIC CHAINS 



11 



singular zeroes between and k + 1 for t G [—So, So]- By Lemma[2J (A3) and then 
(Al) and (A4), we see that for all <C s < So 

z m +i,m+fc(w 1 (e) - u 2 (e)) - z m+ i !m+fe (u 1 (-e) - u 2 (-e)) 

/ 42 n _f k-1 zero of Type I < _k 

^ ' ' \ k zero of Type II _ 2 ' 

Now, given any < e < So, as all the zeroes of w 1 (±e) — u 2 (±e) are regular at 
i = to + 1, m + k, we can by continuity find neighborhoods U 1 and U 2 of u 1 (0) 
and w 2 (0) so that for any w 1 (0),u 2 (0) G U 1 x C/ 2 , 

z m+ i, m+ fc(tt 1 (e) - u 2 (e)) = z m+ i i „ 1+fc (t; 1 (£) - v 2 (e)), 

Zm+i,m+i : (« 1 (-£)-M 2 (-£)) = z m+ i i „ 1+fc (t; 1 (-e) - v 2 (-e)), 

and so that u x (f) — v 2 (£) is not zero on {to, m + k + 1} X [— e, e]. The claim now 
follows from applying Lemma [2j (A3) and then (Al) to v 2 - ti 1 . □ 

The zeroes of infinite degree however do not necessarily persist for small pertur- 
bations in the product topology, as a zero in a small neighborhood can "escape" 
towards ±oo. We can however show the following: 

Lemma 6. Assume u(0),v(0) G /C such that u(0) — v(0) has a zero of degree 
co. Then for each e > there exist |io|,|ii| < £ so that either u(to) — v(ti) or 
u(to) — u(ti), u(to) 7^ u(t\) has a singular zero of a finite degree. 

The proof of Lemma [S] is simple but quite involved as it requires analysis of a 
number of cases with respect to the sign of derivatives du 1 (0)/(ii, du 2 (0)/dt. We 
give it in Appendix 2. 

4.3. The DC case. We now state analogues for the Lemmas above for du/dt 
instead of u 1 — u 2 . Assume in the next two lemmas that F is DC. 

Lemma 7. Assume u(0) G K,, such that du/dt has a singular zero of a finite degree 
k at (m + 1,0),..., (to + k, 0). Then for all sufficiently small e > there exist a 
neighborhood U of u(0) such that for all v(0) G U, 

m-\-k , 

^2 dj(dv/dt;t - s,t +s) > -■ 

j—m+l 

The proof is analogous to the proof of Lemma [5] Similarly, we can show (for the 
proof also see Appendix 2): 

Lemma 8. Assume u(0) G K. such that du/dt has a zero of degree oo. Then for 
each £ > there exist \to\, \ti \ < e so that u(to) — u(ti), u(to) u(ti) has a singular 
zero of a finite degree. 

5. The average number of intersections and zeroes 

This section is the core of the paper. We apply here the results of the previous 
section to solutions of (jl.ip . and define average number of self- intersections and 
derivative zeroes of a translationally invariant measure on X , and average number of 
intersections of two probability measures on X. We then show that these functions 
(denoted by Z) are weak Lyapunov, that is non-decreasing with respect to the 
evolution of (|1.1[) . Furthermore, if there is a singular zero in the support of these 
measures, we show that Z is strictly decreasing. Finally, these properties are in 
some sense continuous in the weak topology of measures. 

These ideas are measure-theoretical equivalents of the properties of the function 
z counting intersections of two solutions of ([l.ip and (|1.2[) ([5]). 
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The intersection-counting function (u,v) i — > Zi,j(u — v) is not invariant for the 
transformation R, thus not well defined on the quotient sets X , /C. To fix that, we 
introduce two-argument functions Zi, Zij : X x X — >N U {0} with 

Zi(u,v) = y~]zj(u - v + r), 

r6Z 

2i,j-(w,u) = -v + r). 

rGZ 

First note that z,, are finite. Indeed, if it, « € /C„, then it is easy to see that 
Zi(u, v) < 2n + 1, thus Zij(u, v) < (2n + — i). Also z,-, 2jj are invariant for the 
transformation R in either argument, thus well defined on X x X. 

Our key object will be the space of T-invariant Borel probability measures A4QC) 
on /C, and 

oo 

M(X) = |J M()C n ). 

k—l 

We will always consider the weak topology on Ai(JC) and M{X) with respect to 
the product topology on X. Note that /C is homeomorphic to a product of compact 
sets, thus compact. As /C is compact and T is continuous, Ai(JC) is closed, and by 
the Tychonoff-Alaoglu theorem A4(K) is compact. 

Before proceeding, let us describe the space ftA{JC). It is clearly non-empty 
(containing e.g. the Dirac measure 8q, where is u = 0), and actually very large. 
Most invariant measures of area-preserving twist maps generated by any generating 
function V are naturally embedded in M(K n ) for sufficiently large n (see Section 
[2]). We can also embed translationally invariant measures on 5 Z , where S is a 
bounded measurable subset of K, in different ways in A4(IC). 

The semiflow (p now naturally induces a continuous semiflow ip\ on M(IC n ). 
The flow (p* on X also induces a local flow ipl on A4(X), as one can easily see that 
one can find S(n) > so that for all u(0) € /C„, u{t) G /C n+ i for t 6 [— <5(n),0] (e.g. 
by compactness argument). We will often use the notation fi{t) — iplfi(0). 

Given fi \fi . ,/i G A^(/C n ), we define the intersection-counting functions 

Z^i 1 ,^) = J J ZQ(u,v)d^ 1 (u)d^ 2 (v), 

Z{p) = J z Q f^j dfj,(u), 

where du/dt is given by the right-hand side of By definition, 

< ^(mSm 2 ) <2n+l, 
< Z(fi) < 1. 

The function Z^ 1 ,^ 2 ) measures the average number of intersections with re- 
spect to two measures; the function Z{p) measures the average number of "self- 
intersections" with respect to and Z(/j) measures the average number of critical 
points duo/dt = with respect to (jl.ip . First note that the definitions above are 
invariant of the spatial location as /z 1 , /i 2 , /i are by definition T-invariant, thus the 
following holds for any integers n < m: 

^(/i 1 ,/! 2 ) = J z n (u,v)d^ 1 (u)dfj 2 (v) 
1 



Zn,m(u,v)dfi 1 (u)dfJ, 2 (v), 



analogously for Z(p), Z(ji). 
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We will see that Z, Z, are "almost" Lyapunov functions in the sense of LaSalle, 
thus they severely constrain the dynamics of ip^,. Generally the properties below 
will hold for functions Z in both AC and DC cases, and for the function Z in the 
DC case only. 

Proposition 2. If/j}(0), p 2 (0), p(0) G M(JC), then the functions t h^- Z(p}(t), p 2 (t)), 
t M> Z(p(t)) are non-increasing. 

Proof. We write shortly p 1 , p 2 instead oi p 1 (0) , p 2 (0) . Let 
(5.1) Z(t):=Z(p\t),p 2 {t)). 
We now have for < s < t 

Z(t)-Z(s) = [ {z (u(t),v(t)) ^ zo(u(s),v(s))}dp 1 (u(0))dp 2 (v(0)). 



Applying Lemma [2 first the zero-balance equation (Al) for m = 0, n = 1, and 
then the T-invariance property (A2) of the crossing functions co,ci and the fact 
that p 1 , p 2 are translationally invariant, we easily get 

(5.2) Z{t) - Z{s) = - ^ J d a (u-v + r;s,t)dp 1 (u(Q))dp 2 (v(Q)) 

which is < 0. The case Z(p(t)) follows by inserting p 1 — p 2 . □ 

The relation (|5.2p will be very useful in the following, as it relates the occurrence 
of singular zeroes with the rate of decrease of Z. For future reference, because of 
translation invariance (Lemma [2J (A2)), for any < s < t, and integers m < n, 
(|5.2|) can be written as (using again shorthand (15. lp ) 



1 f 

(5.3) Z(t) - Z(s) = V V / di(u -v + r; s, tidu^MO^dfStviO)). 

n — m z — ' z — ' J 

i—m r£Z 

We will say that two measures p 1 , p 2 G M.(K.) have a non-transversal intersec- 
tion, if there exist u ,u in their supports so that u 1 — u 2 + r has a singular zero 
for some reZ. A non-transversal intersection is proper, if there is a singular zero 
as above with a finite degree. 

We also say that any non-increasing function t — » f(t) is strictly decreasing at 
to, if for each e > there exists S > so that f(to + e) < /(^o — e) — 6. 

We now show that a proper non-transversal intersection implies that Z(p 1 ,p 2 ) 
is strictly decreasing and that this also holds for small perturbations in the space 
of measures. 

Proposition 3. Assume p 1 (to) , p 2 (to) G M(1C) have a proper non-transversal in- 
tersection. Then there exist neighborhoods U^,\J^ of p 1 (to) , p 2 (to) in A4(IC) so 
that for all ^(to) G U*, v 2 (t a ) G U£, t <-> Z (v 1 (t) , v 2 (t)) is strictly decreasing at 
to- 
Proof. Without loss of generality set to — 0, and choose u 1 , u 2 G K in the supports 
of ^t 1 (0), /^ 2 (0) respectively, so that u 1 — u 2 has a singular zero of a finite degree 
k G N at (m,to), ...,(m + k — l,t ). By Lemma [5J we find open neighborhoods 
t/ 1 , t/ 2 of u 1 ,^ 2 so that (14~T1) holds for any w 1 (0),w 2 (0) G U 1 x U 2 . As u x ,u 2 are in 
the supports of U 1 , U 2 respectively, p}(U x ) > 25 , p 2 (U 2 ) > 2S for some S > 0. 
Let U*, C/ 2 * be the neighborhoods of p 1 (t ) , p 2 (t ) in M(IC) so that for all v 1 G 

(5.4) ^(f/ 1 ) ><5 , v 2 (U 2 ) >6 . 
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Now if Z(t) = Z(v (t), v 2 (t)), combining (I5.3[) (summing over m + 1, m + k, 
with t = e,s= — e) with (|4.ip wc get 

Z(e)-Z( £ ) < -iyi Ir x(t; 1 (0))l cr .(t; a (Q))d/i 1 (« 1 (0))d^ 2 (i; 2 (0)) 

where 1[/ denotes the characteristic function of a set U. □ 

We now extend this to functions Z(p(t)), Z(p,(t)). 

Proposition 4. Assume p-(to) £ M{JC) has a proper non-transversal intersection 
with itself. Then there exists a neighborhood U* of p(to) in M.QC) so that for all 
t/(to) £ U* , t i— > Z(v(t)) is strictly decreasing at to. 

Proof. This follows from Proposition [3] by inserting p, = p 1 = p 2 and taking 

u* = uinu^. □ 

Proposition 5. Assume F in U.l)) is DC. Choose /u(to) € A^(/C) with u £ K, in 
its support such that du/dt has a singular zero of a finite degree. Then t <— > Z(p(t)) 
is strictly decreasing at to . 

Proof. Analogously as Proposition |4j applying Lemma [7] instead of [5] □ 

6. The characteristic flow-map and maps 

This section is dedicated to the proofs of Theorems [1] and [21 Recall that a 
probability measure p on /C is (tp, T)-invariant in the DC case, if it is invariant for 
all (p*, t > (time invariance) and for T (space invariance); and in the AC case 
if invariant for ip 1 and T. We say that a set A C K, is (tp, T)-invariant, if it is in 
the DC case invariant for all y>*, t > and T; and in the AC case invariant for 
ip 1 and T. Keeping in mind this, we will continue to treat the AC and DC cases 
simultaneously. 

As ergodic theory for a commuting semiflow (in the AC case non-autonomous) 
and a homeomorphism is not standard, we give a few comments. First, it is not 
a-priori clear that (ip, T)-invariant measures exist on a compact, (tp, T)-invariant 
set such as JC. We address this in the next section, and construct a rich set of 
invariant measures. 

A (tp, T)-invariant measure \i is ergodic, if for any Borel measurable (tp, T)- 
invariant set A, fi(A) £ {0, 1}. An analogue to the ergodic decomposition theorem 
holds ([H>], Proposition 4.1.12), so we will in some instances without loss of gener- 
ality assume that fi £ M.(K) is ergodic. Note, for example, that given any ergodic 
measure \i £ M(X), by ergodicity there is n large enough so that fi(K n ) = 1, 
thus supp(/it) C IC n as K- n is closed. Therefore assuming /x £ A4(JC) rather than 
H £ M(X) is not a restriction. 

Also note that, given a (tp, T)-ergodic measure, the Birkhoff averages with respect 
to tp* or T do not necessarily p-a.e. converge to the same value. The reason is that 
a (tp, T)-ergodic measure is not necessarily either </>-ergodic or T-ergodic. 

First we show that tp is invertible on the support of any (tp, T)-invariant measure. 

Lemma 9. If fi is (tp,T) -invariant supported on K. n for some n £ N, then in 
the DC case tp t \ suvp ^ is a continuous flow, and in the AC case tp 1 \ supp ^ is a 
homeomorphism. 

Proof. As supp(/i) is (tp, T)-invariant and tp is continuous, it is easy to see that for 
any t > (in the AC case t an integer), tp 1 (supp(^)) =supp(/i). By Lemmad] tp 
is a local flow, thus injective. We conclude that ip* is bijective on supp(/x). As 
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X is metrizable (a quotient space of a metric space, with metric induced by any 
norm (|3.3[) ) and IC n is compact, supp(/i) is compact and ip 1 :supp(/x) — >supp(/it) is 
a homeomorphism for all i. □ 

Proof of Theorem^ Assume /i 6 A4(JC) is (tp, T)-invaxiant. First we show that 
t M> Z(/i(i)) (where fi(0) — /i) is constant. In the DC case, this must be true as 
fj,(t) = fi(0) for all t > 0. In the AC case, this follows from /i(n) = fi(0) for integer 
n, and Lemma [5] which showed that t 1— >■ Z(n(t)) is non-decreasing. 

It is clear from the definition that the map 7r : /C — >T x M defined with (11.41) is 
continuous. We first show that 7r| supp ( M ) is injective. Indeed if it 1 ,?/ 2 esupp(^) are 
such that ^(u 1 ) = ir(u 2 ), then u 1 — u 2 + r has a singular zero of a finite or infinite 
degree for some r £ Z. If the degree is finite, by Proposition [H t i-> Z(v(t)) is 
strictly decreasing at i = 0, where z^(0) = /x, which is a contradiction. Now if the 
degree of the zero is infinite in the DC case, by Lemma [5] and yj-invariance of fj,, 
one can find other v , v 2 Gsupp(/i) so that v — v 2 + r has a singular zero of a finite 
degree, which reduces it to the previous case. 

Consider the case of a singular zero of infinite degree in the AC case. By 
Lemma [HI we find w 1 (to) — v 2 (^i) with a singular zero of a finite degree such 
that w 1 (0), , y 2 (0) Gsupp(//). By the same argument as above, we show that t H> 
Z(n(t + to),[J,(t + ti)) is non-decreasing with period 1, thus constant, which is a 
contradiction with Proposition [3] and v 1 (to) £supp(/i(io)), v 2 {ti) Gsupp(/j,(£i)). 

Now assume the contrary, i.e. that there are two (ip, T)-invariant measures 
(Jr,fJ, G A4(IC), so that for some u l Gsupp(/i 1 ), u 2 esupp(/i 2 ), we have ^(u 1 ) — 
ir(u 2 ). But [i — /i 1 1 '2 + /i 2 /2 is also a (<p, T)-invariant measure whose support 
contains both u , u 2 , which is impossible. 

Finally, the case /i £ A4(X) reduces by the ergodic decomposition theorem to 
analysis on M(lC n ) for sufficiently large n. □ 

Remark 1. We actually proved more in the proof above: if u,v £ A, then u — v 
can not have a singular zero. 

Proof of Theorem^ By definition, it : JC — >• /C is continuous and by Theorem[T]it is 
bijective. As K, is compact and /C a metric space, 7r : /C — » K, is a homeomorphism. 
The claim now follows from Lemma [5] and definition p.5p of T, (p. □ 

We complete the Section with a conjecture and a couple of comments on the 
likely structure of the characteristic map T associated to (11.11) . 

Conjecture 1. The map T can be canonically extended to a homeomorphism T : 
T 1 x R 4 T 1 x 1, diffeomorphic to an area-preserving twist diffeomorphism fv,F- 

The numerical simulations [7] and our construction suggest the following descrip- 
tion of T. Given an AC forcing F, we can define 

F = [ F(t)dt, 
Jo 

a(F) = fj\F(t)-F) 2 dt 

We conjecture that for a fixed V,a(F) (including the DC case cr(F) = 0), there 
is F c such that for all F > F c , the map fv,F is integrable. Conversely, for a fixed 
V, F, the limit cr(F) oo corresponds to an " anti-integrable limit" of the family 
of maps fv,F- 




16 



SINISA SLIJEPCEVIC 



7. The Aubry-Mather theorem 

This section is dedicated to the proof of Theorem [31 We first define rotation- 
ally ordered configurations in R z and the rotation number and deduce their main 
properties, then show in general existence of (tp, T)-invariant measures, and from 
this complete the proof. 

Recall the definitions of translations T(u)i = ztj+i and Ru = u + 1 commuting 
with (jl.ip and each other. We can define T pq = T p R q , or T P}q (u)i = Ui +p + q for 
any p, q G Z. We say that a configuration u G M, z is rotationally ordered, if for any 
p, q G Z, either T p ^ q u < it or T pq u > u (i.e. the set of all T pq translations is totally 
ordered). 

The following is an elementary result of Aubry-Mather theory; for a proof see 
e.g. Bangert [BJ, Corollary 3.2: 

Lemma 10. If u G M. z is a rotationally ordered configuration, then there exists a 
unique p G K so that for all i,j G Z, 

(7.1) |iij -tij - p(j - i)\ < 1. 

The number/? can be characterized as p = sup{p/g : T p ^ q u < u} = ini{p/q : 
T Ptq u > u}. 

If u G M. z is such that the limit 

p{u) = lim 



|i-i|->oo j — i 

exists, we call it the rotation number (or mean spacing) of a configuration. For 
rotationally ordered configurations, pin (|7.1I) is its rotation number. 
Now we show existence of (ip, T)-invariant measures. 

Lemma 11. Assume 1Z C K, is a closed, non-empty (<p,T) -invariant set. Then 
there exists a (ip,T) -invariant probability measure supported on 1Z. 

Proof. Choose any Borel probability measure v on 1Z (say Dirac supported on any 
element). We in the DC case define 



\ ' 



<piT*v dt, 



and in the AC case 



3=1 fe=l 

Here p\ , is pulling of measures with respect to p t ,T k . As 7Z is compact and v n 
are well defined probability measures supported on 1Z, v n has a weakly convergent 
subsequence v nk —t [i and fi is supported on 1Z. It is easy to check that /i is (ip, T)- 
invariant (e.g. by checking this for continuos functions / : 1Z — KR and then applying 
the Riesz representation Theorem). □ 

Proof of Theorem [3J Denote by 1Z P the set of all rotationally ordered configurations 
with the rotation number p. It is non-empty, as u defined with Ui = pi is in 
IZp. From (|7.1[) one immediately sees that 7Z P C IC n where n is any integer > 
\p\ + 1. From the definition it follows that 1Z P is closed, thus it is compact. As T p ^ q 
commutes with ip and T, by Lemma |3]<y9 is order preserving and T is trivially order 
preserving, we see that 1Z P is (p, T)-invariant. It is now enough to apply Lemma 
HUto TZ P . □ 
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8. Description of invariant measures in the DC case 

Here we prove Theorem^! and so describe the set of (tp, T)-invariant measures on 
X and /C, denoted by £{X) and £ (K.), in the DC case. The standing assumption in 
this section is that F in (jl.ip is DC. We first define the average speed of a measure 
and a configuration. We then show by applying the Poincare-Bendixson theorem 
that any /i E £{K) is entirely supported on equilibria and periodic orbits. If the 
average speed is ^ 0, then \x is supported on a single periodic orbits. Finally, if the 
average speed is zero, we show that /x is supported on equilibria. 

Let /j, £ £{X) and u £ X. We can define their average speeds v(fi), v(u) as 



(8.1) v( f i) = I (uo(l)-«o(0))^, 

/o o\ ,„ \ r u o(t) ~ «o(0) 

(8.2) v(u) — hm , 

t— s-oo t 

where v(u) is defined when (|8.2j) is convergent. 

Lemma 12. //// is (ip,T)-ergodic on X, then for /i-a.e. u, v(u) — v(p). 

Proof. As /x is ^-invariant (but not necessarily i/3-ergodic!), there is a set of full 
measure such that v(u) is well defined. As by definition and Lemma 0] the spacing 
\ui + i(t) — Ui(t)\ is bounded uniformly in i,t, the definition (|8.2[) is independent 
of i in Uj, thus v(u) = v(Tu). Thus v(u) is defined on a set of full measure, tp 
and T-invariant measurable function, so it must be constant /x-a.e. and equal to 
v(p). ' □ 

We will use a version of the well-known Poincare-Bendixson theorem on asymp- 
totics of 2-dimensional flows. For a general discussion of the Poincare-Bendixson 
theorem see e.g. [TB], Section 14.1; for the proof of the exact formulation below 
see [9], the proof of Proposition 1 (the only difference being that they consider 
compact subsets of the plane, which makes no difference, as the argument works 
for any compact sets which can be embedded in a two dimensional sphere). 

Lemma 13. Assume K is a compact subset of the cylinder T 1 xl, and <p a complete 
continuous flow on K generated by a C° vector field on K. Then the only recurrent 
points of ip are equilibria and periodic orbits. 

We need the next lemma to be able to claim that if tt(u) is an equilibrium, then 
u must be an equilibrium. 

Lemma 14. If u £ A, then w = du/dt can not have a singular zero. 

Proof. Find a measure fi £ £()C) for some K, — /C„ so that u is in the support of /j,. 
As /i = /i(0) is by definition ((p, T)-invariant, t i-> Z(n(t)) is a constant function. 
Now if w = du/dt has a singular zero of a finite degree, from Proposition[S]we obtain 
that Z(fj,(t)) is strictly decreasing at t = 0, which is a contradiction. If w has a 
singular zero of infinite degree, by Lemma |S] we can find v 1 , v 2 in the support of /i 
so that w 1 — v 2 has a singular zero of finite degree. This implies that by Proposition 
[3l t n> Z(fi(t)) is strictly decreasing at t = 0, which is again a contradiction with /i 
being (ip, T)-invariant. □ 

We now show that all of this combined with our results on the projection 7r gives 
a fairly good description of supports of \x £ £(X). 

Lemma 15. Assume fi is ((p,T)-ergodic. Then /x is supported on ^-periodic orbits 
with the same period to > 0, or on equilibria. 
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Proof. As [i is (ip, T)-ergodic, there exists K, — JC n for n large enough so that 
fi G £(JC). Take K, (p. as in (|1.5p . Now we deduce from Theorems [T] and [3] that 
most of the assumptions of Lemma [TBI are satisfied; we only need to show that ip is 
generated by a C° vector field. One can easily check that for (x,p) G T 1 x R, the 
vector field is given by 



(8.3) 



(x,p) l-> 



d(ir 1 (x 1 p))o/dt 

d(ir~ l (x , p)) \ / dt — d(ir~ 1 (x,p))o/dt 



where duo/dt, du\/dt for u = tt^ 1 (x,p) are defined with the right-hand side of 
(II. ip . From Lemma 1131 we deduce that all ^-recurrent orbits on K are periodic 
orbits and equilibria. Consider the pulled ^-invariant measure 7r*/i on K. As K 
is metrizable and compact, the set of ^-recurrent points has full measure with 
respect to n^fi (see e.g. [TB], the proof of Proposition 4.1.18). Now note that by 
Theorem [I] ir^ 1 maps periodic orbits into periodic orbits, and by Lemma |141 ir^ 1 
maps equilibria into equilibria (otherwise the right-hand side of (|8.3[) would be zero, 
i.e. for u = tv (x,p), du/dt would have a singular zero, which is impossible by 
Lemma [T4l . If we denote by Vt the orbits in /C with period t, t > (where Vo are 
equilibria), we now have 



(8.4) 



v t>o 



Now for each t >0,Vt is (<p, T)-invariant, thus by ergodicity ^(Ut<t 'Pt) G {0, 1}. 
As to i — > /i(Ut<t 7- , t) is non-decreasing, there exists a unique to so that p>(Vt ) = 1. 
As Vto is closed, it contains the support of /i. □ 



We finally need to show that, if \i is supported on periodic orbits and has non-zero 
speed, it is supported on a single one. 

We say that the orbit of u is totally ordered, if the set T(u) = {</? t T fc u, t G M, k G 
Z} is totally ordered. That means that for any v,w G T(w), we have either v > w 
or v < w. 

Lemma 16. Assume /i is [tp,T)-ergodic, v(p) ^ 0, and supported on periodic orbits 
with period to > 0. Then it is supported on a single totally ordered periodic orbit 
such that the tt o T J -image of supp(fi) is independent of j G Z. 

Proof. We prove it in three steps. First we show that tt(io) = u(0) ± 1, then we 
show that u(0) and u(t) in the support of jj, can not intersect, and finally we show 
uniqueness of the periodic orbit. 

The average non-zero speed is by Lemma [T^] the same for all the periodic orbits 
in supp(/i). By Theorem [I] it projects each of them to a homotopically non-trivial 
invariant circle on T 1 xK (otherwise the speed would be zero), thus u(to) = u(0)±l, 
with the sign the same as the sign of v(fi). Assume without loss of generality 
u(t ) = u(0) + 1. 

Now assume u(0) and u(ti) intersect, < t-y < to (the only analogous alternative 
is — to < ti < 0). Consider w(t) — u(t) — u(0). Recall that by Theorem [1] and 
Remark[TJ all intersections in the support of [i are transversal, so w(t) has a regular 
zero, say at i. Then w(t) intersects zero at 

(8.5) x(t) = + 

Ui - U l+ i 
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As w(t) has no singular zeroes, it is easy to see we can continue the function 
t M> x(t) locally with (|8.5p if x(t) not an integer, and if x(t) = i an integer, with 



x(t) 



Thus we can find a continuous function 1 1-> x(t) so that w(t) has a regular zero 
at L^WJ ■ But by assumptions w(to) = 1, so the only possibility is that for some t 2 , 
t± < t 2 < to, lim^j- a;(i) = +oo or lim t t - = — oo. Assume without loss of 
generality lim^j- x(t) = +oo. As by the proof of Lemma[|J t h-> w(i) is continuous 
in /qo, we easily see that lim^oo Wi(t 2 ) = 0, thus 

(8.6) lim \ Ui (t 2 ) -Uj(0)| = 0. 

By compactness, we can find a sequence T™ fc , — »■ oo such that T 7lk u(0) converges 
to some v. As T commutes with tp, from (|8.6|) we see that = w (0) ; and f is in 
the support of ft. But t2 < to, which is in contradiction with the assumption. 

Now we see that each periodic orbit in the 7r-image of support of fi intersects 
each x coordinate exactly once at a point p, where (x,p) e T 1 x 1. We can thus 
identify the quotient set supp(/i)/V with a compact set B of intersections at x = 0, 
and consider B to be a subset of R parametrized with p. By definitions, T induces 
a homeomorphism T on B, and the induced flow tp is constant on B. The induced 
measure Jl is (T, (^)-ergodic on B, and as ip is constant, it must be T ergodic. Finally, 
note that T is order preserving on B, that means that p > q implies T(p) > T(q). 
We see that, as otherwise we find u,v in equivalence classes of p,q which violate 
order preserving, such that Uo(t) = Vo(t'), and then u(t) — v(t') has a singular zero 
at 0. As by assumptions, B is the support of p,, T a homeomorphism and p. is 
T-ergodic, there must be po £ B such that its w-limit set with respect to T is the 
entire B. But T n (po) is monotone and bounded, thus converges to a single point. 
We conclude that B must be a single point, so support of [i consists of a single 
periodic orbit. 

It must be a totally ordered set, as we have shown that u(i), u(t') do not intersect 
unless equal. □ 

Proof of Theorem^ Take any (tp, T)-ergodic measure fi on X. If = 0, con- 
sider the function J V(u , u{)dfi. The derivatives of V and u, are by compactness 
uniformly bounded on the semiorbit of u(t). Thus we can exchange the integral 
and derivation, use T-invariance of fj,, (jl.ip and finally (p-invariance of fx and obtain 

^ J V(u ,ui)dfi = J (Vi(u ,ui)u' (t) +V 2 (uo,ui)u' 1 (t))d(i 

(Vi(u Q ,Ui)u' (t) + V 2 (u- ll u )u' (t)) dfi 
(u' (t)f dfi + F J u (t)dfi 
(u' (t)f dfi + F J Qf u' (t)dtj dfi 

(8.7) = -J (u'^fdfi + Fvtfi). 

However V is continuous on X and fi is ^-invariant, so the left-hand side of (I8.7[) 
must be 0. We deduce that for /i-a.e. u, u' Q (t) = 0, so by continuity argument and 
T-invariance we see that the support of fi consists of equilibria. 

If v(fi) ^ 0, the claim follows from lemmas [T5l and [T6l □ 
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Remark 2. The relation 



F-v(fi) = / (u' (t)) 2 dp 



obtained above, thus valid for DC dynamics, has standard physical interpretations, 
such as forcexspeed= energy dissipation, or voltagex current—resistance. 



In this section we further discuss the DC case when the average speed is not 
zero. We show that the periodic orbit is then a uniformly sliding state, and discuss 
the properties of its " modulation function" . The results of this section have already 
been obtained by Baesens, MacKay and Qin ([3], [5], [32], [33]), by combination of 
a limiting procedure (typically going from rational to irrational rotation numbers) 
and a fixed point argument. 

We show that these results follow directly from our ergodic-theoretical and topo- 
logical considerations, hoping for extension of these results to AC driving and sys- 
tems similar to (jl.ll) . 

We say that a solution of is a uniformly sliding solution with the rotation 
number p, speed v and time shift a, if there exists a modulation function m : R — > R, 
m(x + 1) = m(x), such that for each j G Z, t G R, 



Corollary 1. Assume that F is DC, p is (ip,T)-ergodic and v(p) ^ 0. Then p is 
supported on a single uniformly sliding solution. 

Proof. Assume without loss of generality v(p) > 0. As by Lemma \W\ and Theorems 
[TJ [3J 7r projects Uj(t) to a single periodic orbit with period to independently of j 
and with the same projected flow (p, there exists a function rh : R — )• R and a.j G R 
such that fhit + to) = rh(t) + 1 and for all j € Z, Uj(t) = m(t + ctj). It is easy to 
see that v = w(/i) = l/t . We set m(t) := rh(t ■ to) — t, and thus m(t + 1) = m(t) 
and 

(9.1) u j(t) = m{vt + voij) + vt + va.j. 

As the set {u(t),t G R} is by Theorem 2] totally ordered, u(t) is rotationally 
ordered, so by Lemma [TU1 (|7.1|) must hold with some p = p(u(t)). As t i — > 
u(t) — u(0) is continuous in ^(Z), all u(t) must have the same rotation number 
equal to the rotation number of the measure defined as 



where in the second row we used T and y>-invariance of p,. By Lemma fTBl J Q (uj+i (t) — 
Ui(t)) is independent of i as it is the second coordinate of the projection 7r inte- 
grated over one period, and by periodicity independent of u G supp(/i), thus for all 

i G Z, t G R, 



9. Uniformly sliding states in the DC case 



Uj (t) — j p + vt + a + m(J p + vt + a). 




(9.2) 




Combining (|9.1[) and (|9.2I) we get vctj = jp + a for some a G R, which completes 
the proof. □ 
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From Theorems [3J 2] now follow the results from [?] , [32] that if there is no 
equilibrium with a given rotation number pei, there must be an uniformly sliding 
solution with the same rotation number. Furthermore, one can prove continuous 
interdependence of v, p, F and other parameters of (jl.ip by considering continuity 
of (ip, T)-invariant measures and Theorems |31 SI 

10. THE SPATIO-TEMPORAL ATTRACTOR 

In this section we prove Theorem[SJ We need to show that, given any T-invariant 
probability measure ponX and any open set neighbourhood U of X, (|1.6[) holds. As 
lC n are closed and T-invariant, by the ergodic decomposition theorem it is sufficient 
to prove that for fx G M.(fC n ) (denoted again with M.{tC)). 

Recall the definition of P M ,s in Section 2. It is easy to see that the left-hand side 
of (|1.6p can be written as 

(10.1) P„,s(i(u, t), <p*(u) G A}) = I j\^{A))dt. 

Denote the probability measure on the right-hand side of (|10.ip with v s . It is easy 
to see that v s G A4(/C). We want to show that, given any open neighbourhood U of 
A, then lims^oo v s (U) = 1. Assume the contrary, and find a sequence Sk — > oo and 
S > so that u Sk (U) < 1 — 5 for all k. As M(JC) is compact, v Sk has a subsequence 
for simplicity denoted also with v Sk , which converges to some /i G M (K.) . It is 
easy to deduce from (llO.ip that fi is also (^-invariant. Then fi must be supported 
on A, thus fi(A) =1, fi(U) = 1 and as U open, lim^oo v Sk {U) = 1 which is a 
contradiction. 

11. Appendix 1: on the zero-set of a system of infinitely many 

cooperative equations 

11.1. The proof of Lemma [TJ Recall the setting and definitions from Section |3j 
and assume in the following all the assumptions of Lemma [T] 

Lemma 17. Assume w\{§), iffe(O) = is a singular zero of degree k. Then there 
exist real numbers d\, ...,dk such that for all j — 1, k, 

(11.1) Wj(t) = djt j ' + o(t j *), 
where j* = min{j, k + 1 — j}. 

Proof. We prove inductively in I, 1 < I < (k + l)/2, that for all j, I < j < k + 1 — I, 

(11.2) Wj(t) =dj,r t l +o(t l ). 

For I = 1 and all 1 < j < k, as Wj is C 1 and Wj(0) — 0, by the Taylor formula, 

, . dwj(0) , . 

thus 

(11.3) d j ,i = dw j (0)/dt 

We now assume (jll.2l) holds for some I — 1. Inserting (jll.2l) in the right-hand 
side of (|3.1I) and using p.2p (continuity of aj, bj, Cj) we get 

j t Wj{t) = (aj(t)^_i,i_i + 6 i (t)d,-+i >1 _ 1 + ^(tJdjM-O • + ait 1 " 1 ) 

= {aMdo-U-i + bj(0)dj+i,l-i + c^d^t 1 - 1 + o{l)t 1 - 1 + o{t 1 - 1 ) 

(11.4) = (o 3 -(0)d,-_i,,_i + 6^(0)^+1,1-1 + ^(Ojdj.i-i)* 1 - 1 + 
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Integrating it and using Wj(0) = we obtain (|11.2p for I and complete the 
inductive step with 

(11.5) d Jtl = y (^(0)^-1,1-1 + &j(0)d i+ i,,_i + Cj(0)d A i_i). 

□ 

We now need to show that dj ^ and determine their sign. Denote by Sgn(ar) G 
{-1,0, 1} the sign of a; G M. 

Lemma 18. Assume all as in Lemma \T 7 ^ Then for all j = l,...,k, except in the 
case j — k/2, k even, the zero of Type II, we have 

(11.6) sgn(dj) = sgn(w b (0)), 
where b G {0, k + 1} such that \j - b\ < \j - (k + 1 - b)\. 
Proof. From (fTO]) and iJOJ we get 

d M = dtui(0)/dt = oi(0)io (0), 
dfc,i = dw k (0)/dt = b k (0)w k+1 (0), 
dj : \ = for j = 2, k — 1. 

Iteratively from (jll.5[) we easily get using again j* = min{j, k + 1 — j} 

dj = MO) IlLi a<(0)/i*l for j < fe/2, 

^ = Wfe+1 (0) fit, &i(0)/i*! for j > k/2, 

dj = Mo) IIU Oi(o)/i*! + ^+i(o) n- =J fti(o)/i*i for j = fe/2. 

which completes the proof. □ 

Relations (111.11) and (jll.6p now complete the description of the behavior of a so- 
lution around a singular zero of degree k for small t. We summarize it, again assum- 
ing without loss of generality (i.e. up to reparametrization) that (1, 0), (k, 0) G S 
is a zero of degree k at t — 0. In that case, lemmas IT71 and IT51 imply that, assuming 
all as in Lemma [T71 that there exists 5q > such that for all < 6 < So, we have 
(see Figure 1): 

Type of zero k z ,fc+i(w(-<5)) z ,fc+i(w(0)) z 0>k+1 (w(6)) 
I even k + 1 k 1 

I odd k k 1 

II even k k 
II odd fc + 1 fc 



Figure 1: Types of singular zeroes of finite degree 
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We conclude that the number of zeroes is strictly decreasing at singular zeroes. 
By continuity, the number of zeroes does not change at regular zeroes, which com- 
pletes the proof of Lemma [TJ 

11.2. Zeroes of degree oo. As noted earlier, zeroes of degree oo require special 
attention. 

Lemma 19. Assume all as in Lemma]]] Suppose w(0) has a zero of degree oo, 
wq(0) j^z and Wi(0) = for all i > 1. Then there exist dj having the same sign as 
Wq(0) such that for all j > 1, 

(11.7) Wj(t) = djt j +o(t j ). 

Proof. An analogous analysis as in the proofs of Lemmas [T7] and [15] easily yields 
([TT77|) with 

dj = — — J_J_Oi(0). 

■'' i=l 

□ 

11.3. Proof of Lemma [2j We construct explicitly the functions Ci(s,t),di(s,t). 
We first define sets of zeroes S + , S~ which are subsets of the set of singular zeroes 
S C Z x [Tb,Ti] which "disappear" either at t = to or at t — t^. Assume without 
loss of generality that (1, to), (k, to) is a singular zero of degree k. Then we define: 

(i) Zero of Type I, k even, zero at i = 1, 21 = k: 

S~ 3 (Mo), 

S+ 3 (l,to),(2,to),..,{l-l,to),(l + l,to),...,(2l,t ). 

(ii) Zero of Type I, k odd, zero at i = 1, 21 + 1 = k: 

S+ 3 (l,io), (2, to), (I, to), (I + 2, to), (21 + Mo). 

(iii) Zero of Type II, k even, zero at i = 1, 21 = k: 

S+ 3 (l,t ),...,(2l,t ). 

(iv) Zero of Type II, k odd, zero at i = 1, 21 + 1 = k: 

S- 3 (l + l,t ), 

S + 3 (Mo),-, (2/ + Mo). 

(v) Zero of degree oo at i > 1 : 

S+ 3 (1, t ), (2, to),- 

Also S~ , S + contain no other elements apart from those defined above. We now 
define 

di(s,t) = \S + n {i} x [s,t)\ + \S~ n {i} x (s,t]|, 

where \A\ is the cardinal number of a set. The rationale for this is clear from Figure 
1 and ([TTTT)l . (dTTJ). 

Without loss of generality we construct co(0,S) for some S > 0. For each 
t £ [0, S 1 ] we choose an open interval J t Cl containing t in the following way: 

(i) If wo(t) ^ 0, then I t is such that for all s £ I t , wo{s) ^ 0. 

(ii) If wo{t) — 0, a regular zero, then 7 t is such that for all s £ I t , Wi-i(s) ^ 0, 
m+i (s) 7^ 0. 

(iii) If wo(t) = is a singular zero of Type I, degree odd, zero at — /, — I + 1, / 
(i.e centered at 0), we choose I t so that for all s £ Jj, s^f, we have %(s) ^ for 

i = -i,-, -1,1,-, I- 
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(iv) If If Wo(t) = is a singular zero in all other cases, It is such that for all 
s £ It, s t, we have wq(s) ^ 0. 

We constructed the neighborhoods It in such a way that for each time s E It, 
there is at most one zero which may be crossing i = (the singular zero cases follow 
from all the discussion above in this section.) 

We first define c(s,s') for w(s), uu(s'), s < s' in the same interval It, by cases 
(i)-(iv) above: 

Case (i) co(s, s') — 0, 

Case (ii) c (s, s') = z (w(s')) - z (w(s)), 

zo,i{w(s')) + 1 - zo,i(w(s)) s'>t>s 
otherwise, 
Case (iv) co(s, s') = 0. 

Clearly cq(s, s') E { — 1, 0, 1}. Now assume It ,It! , ■■■It j is a finite subcover of [0, S], 
and = so < si < ... < s^. = S such that s,_i, Si are in the same open interval I ti . 
We define 

r-i 

c Q (0,S) = ^2 c (s l -i,s i ). 

It is now straightforward to verify that the definition is independent of the choice 
of ti, and that Ci,di satisfy (A1)-(A4). 

11.4. Proof of Lemma [3j We slightly adapt 5 , Proposition 2.1. Let C be the 
uniform bound on |cj(f)|, and let Wi(t) — exp(Ct)wi(t). Then Wi(t) is a solution 
of an equation of type (|3.ip with 5j = aj,6j = 6i,c,; = Cj + C, thus ai,bi,Ci > 
satisfying p.2p . Write it shortly 

(11.8) <Mv/dt=A(t)w(t), 

where A(t) is an infinite matrix. By standard results on existence of solutions on 
Banach spaces, the solution of (|11.8[) exists locally in any of the norms (|3.3I) and is 
the limit of Picard iterations 

w^(t) = w(0)+ f A(s)w^ 1) (s)ds, 
Jo 

vj(°\t) — W(0). As A(s) has non-negative elements, we easily see that w(0) > 
implies w(t) > for t > 0, and Wi(t) and Wi(t) have the same sign. 

12. Appendix 2: On zeroes of degree oo 

This Appendix is dedicated to proofs of Lemmas [5] and |H] Assume without loss of 
generality that the zero of degree oo of w(0) = u(0) — v(0) is at i = 1, i.e. u>o(0) ^ 0, 
Wi(Q) = for alH > 1. We first consider the case when all Ui are strictly monotone 
on either sides of t = 0. The proof is somewhat combinatorial: we consider different 
alternatives for the signs of duj/dt. 

Lemma 20. Assume all as in Lemma\(% Assume also that for some iq > 1, there 
exist Sj > so that duj/dt is not zero on ( — <5j- , 0) and {0,Sj) for all j > iq. Then 
the conclusion of Lemma\^ holds. 

Proof. First note that the conclusion of Lemma [T^l holds for w(0) = u(0) — v(Q), 
and recall that %(0) = Vj(0) for j > 1. 

Case 1. Assume that for some ji > iq, the signs of dv,j 1 /dt,...,duj 1 +3/dt are 
+,+,+,+ (analogously we show — ) forte (— 5j l: 0), (— 8^+3, 0). Now 

because of (|11.7[> . for some e > small enough, for any to G (—£,0), Uj x (to) > 
v n (t ), u jl+1 (t ) < v jl+ i(t ), u n+2 (t Q ) > v jl+2 (t Q ); or u jl+1 (t ) > v jl+1 (t ), 
Uj 1 +2(to) < i>j 1 +2(io), u ji+3{to) > v ji+3(to)- Consider without loss of generality the 



Case (iii) cq(s, s') = 
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first possibility and fix some > to > —e. Now find > t\ > to so that = 
Vj 1+ i(to). As by assumptions Uj 1 (ti) > v^fo) and Uj 1 +2(^1) > ^1+2(^0), we see 
that u(ti) — v(to) has a singular zero of degree 1 at ji + 1. 

Case 2. Assume that for some io < ji < 32 < 33 < 3 a (not necessarily 
consecutive), the signs of duj 1 Jdt, ...,duj 4 /dt alternate on (0, 5j x ), (0, <5j 4 ), say 
— ,+,—,+ (analogously —,+,—,+). By (|11.7p . either Uj(t) > Vj(t) for all j > io 
and all t £ (0,<5j), or < i>j(£), consider again only the first possibility. 

Fix < to < s for e small enough, and then we can find < <2 < to so that 
^2(^2) = Vj 2 (to). Thus we can find t% < t\ < to, as the smallest such time such 
that Uj(t\) > Vj(to) for all j = j\ + 1, j'3 — 1, and then for t — t\ at least one of 
the inequalities is =. As by assumptions Uj x {t\) > Vji(to)> u j 3 (ti) > v j 3 {ta)i we see 
that u(ti) — v(to) has a singular zero of finite degree. 

Case 3. Assume that for some j\ > io, k > 2 the signs of duj 1 /dt, ...yduj^k/dt 
have the signs — , +, +, +, — on (— Sj, 0) and the signs — , — , — on (0, <5,-) (or 
the alternate sign case). Fix < to < £ for e small enough, and find — e < t\ < as 
the largest t < so that Uj (t) > itj (to ) for all j = j\ + 1, ... , j± + k — 1, and then for e 
small enough at least one of inequalities is =. As by assumptions, Uj x (ti) > Uj 1 (to), 
u ji+k{t\) > Uj 1+ k{to), u(ti) — u(to) has a singular zero of a finite degree. 

Now it is straightforward to check that one of the cases 1-3 (including the alter- 
nate possibilities) must hold. □ 

Proof of Lemma O We consider separately the DC and AC case. In the DC 
case, we use the fact that all the results on the zero-set apply to w = du/dt. As 
this does not necessarily hold in the AC case, we apply real analyticity of the 
solution. 

Consider the DC case. If duj(0)/dt ^ 0, then Uj is strongly monotone in a 
neighborhood of t — 0. Also if duj(0)/dt = 0, by Lemmas [T71 [T51 and [T51 duj/dt — 
djt k + o(t k ) for some dj ^ and an integer k > 1 in all cases except possibly 
when (j, 0) is a part of a zero of Type II of an even degree of w = du(Q)/dt. If 
w = du(0)/dt has only finitely many zeroes of Type II, even degree, on (i,0) for 
i > 1, we see from this discussion that for some i > io the assumptions of Lcmma [20l 
hold and the proof is completed. Now assume the contrary, and consider a zero of 
Type II, even degree of w = du(0)/dt, say at (j + 1, 0), (j + k, 0). Then duj(0)/dt 
and duj + k+i{0) /c?t have alternate, non-zero signs. As there are infinitely many such 
zeroes, we take any three and obtain the Case 2 in the proof of Lemma I2TH 

Consider now the AC case, and assume without loss of generality F(t) in (jl.ip 
is not constant (if it is, this reduces to the DC case). By assumptions and Lemma 
21 Uj(t) is real analytic on every open subset of the set of times for which u(t) £ X. 
Then we have cither 

(12.1) duj/dt = djt k + o{t k ), k > 

or duj/dt — in a neighborhood of t = (see e.g. |17j . Remark 1.2.13). If for 
a ll j > io, the first possibility holds, we can apply Lemma 1201 and complete the 
proof. Assume the contrary, i.e. existence of infinitely many j > 1 so that duj/dt 
locally vanishes around t = 0, and discuss all the possibilities. 1) If there are 
three or more consecutive vanishing duj/dt, inserting it into the derivative of 
implies that fit) is locally constant, thus as real analytic and periodic, constant 
everywhere, which is a contradiction. 2) Assume there are exactly two consecutive 
vanishing diij+i/dt, duj+2/dt, for some t £ (—e,e), where e > is small enough 
such that diij/dt,duj+3/dt have no zeroes on (0,e) (this is possible as they satisfy 
(|12.1[) ). Then for any < to < ti < e, clearly u(to) — u(ti) has a zero of degree 2 
at j + 1, j + 2. 3) Assume there is a vanishing duj + i/dt on t £ (— e, e) such that 
duj/dt and duj + 2/dt have no zeroes for t £ (0,e) and have the same sign. Then 
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analogously for any < to < t\ < e, u(t ) — u(ti) has a singular zero of degree 1 
at j + 1. 4) Finally, consider infinitely many vanishing duj/dt so that the signs of 
duj-i/dt, diij+i/dt are not zero and alternate for t € (0, Sj). But this is the case 2 
from the proof of Lemma l20l As we have exhausted all the possibilities, the proof 
is completed. □ 

Proof of Lemma\^ We can apply Lemma fT9l to w — du/dt and immediately see 
that the Case 3 of the proof of Lemma I2TJ1 holds. So the conclusion of the case 3 
must hold, which completes the proof. □ 
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